Models of Superconductivity in Sr 2 Ru04 
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Recent experimental data on purest S^RuCU single crystals clearly indicate the presence of nodes 
in the superconducting order parameter. Here, we consider one special p-wave order parameter 
symmetry and two two-dimensional /-wave order parameter symmetries having nodes within the 
RuC>2 plane. These states reasonably describe both specific heat and penetration depth data. We 
calculate the thermal conductivity tensor for these three states and compare the results with recent 
thermal conductivity data. This allows us to single out one of these states being consistent with both 
thermodynamic and thermal conductivity data: the planar /-wave state having Bi 9 x E u symmetry. 



PACS: 74.70.Dd, 74.20.-z, 74.25.Bt, 74.25.Fy 
I. INTRODUCTION 

The recently discovered superconductivity in Sr2Ru04 
has been interpreted in terms of a p-wava.jtriplet super- 
conducting state having a full energy gapEla. For exam- 
ple the spontaneous spin polarization seen by muon spin 
rotation experimental and the flat -Knight-shift seen by 
nuclear magnetic resonance (NMR)ta are consistent with 
spin triplet pairing. However, recent specific heat datacl 
and the superfluid densityQ of purest single crystals of 
Sr2Ru04 with T c < 1.5K clearly show low temperature 
behavior consistent with nodes in the order parameter 
very similar to observations of <i-wave superconductivity 
in the high-T c cuprate superconductorsQtl. 

Here we shall study three examples of two-dimensional 
(2D) superconducting order parameters with spin triplet 
pairing having nodes within the RuC>2 a-b plane. The 
first one is the anisotropic p-wave state proposed by 
Miyake and Narikiyotl with A(/c) oc sin^a) ± i s'm(k y a). 
Here, a denotes the lattice constant of the RuC>2 square 
lattice. In order to have a node with this state, how- 
ever, we have to stretch the Fermi wavevector kp towards 
the particular value of kpa = 7r, while a more realistic 
value would-be kpa = 2.7 as judged from bandstructure 
calculations^! In the following we will denote this par- 
ticular p-wave state as the nodal p-wave state. As the 
second and third example we consider the planar /-wave 
states recently proposed by Hasegawa et alO. Here, the 
angular <f> dependence of the order parameter is given 
by A(k) oc cos(2(?!>)e ± ^ and A(k) oc sin(20)e ±l *, respec- 
tively. 

Within circular symmetric weak-coupling BCS theory 
one immediately realizes that the thermodynamics of 
the latter two states is identical to the one of d-wave 
superconductors!! We have worked out the thermody- 
namics of the anisotropic, nodal p-wave state here as well. 
In Figs. ^ and we show our results for the temperature 
dependence of the specific heat C s /jT and the superfluid 
density p s (T) for the nodal p-wave and the 2D /-wave 
states together with the experimental data. For com- 
parison, we also show the results of a 3D /-wave state, 



considered by two of us recently£3. As is readily seen, 
the 2D /-wave states give a better description of the ex- 
perimental data than the 3D /-wave or the nodal p-wave 
states, though the differences between the 2D /-wave and 
the 3D /-wave states are rather small. 




FIG. I. The specific heat C a /-yT as a function of T/T c for 
the 2D /-wave states (solid line) and the nodal p-wave state 
(dashed line) considered in this workpi Also shown are the 
experimental data by Nishizaki et al.u (circles) and the 3D 
/-wave state considered in Ref. |f2| (dotted line). 

Very recently the thermal conductivito-of Sr 2 Ru04 in a 
planar magnetic field has been studiedOllil. Both groups 
studied the thermal conductivity parallel to the a-axis in 
a magnetic field within the a-b plane in a direction tilted 
by an angle 9 from the heat current. Both groups found 
no appreciable angular dependence. This experimental 
result is already inconsistent with the isotropic p-wave 
statc_having a full enery gap and the 3D /-wave state 
statcEj. Indeed, we shall show that the thermal conduc- 
tivity data is consistent with only one of the three nodal 
states considered here: the 2D /-wave state with angular 
dependence cos(20)e ±l< ^. The two other states exhibit 



1 



rather large angular dependence and therefore are incon- 
sistent with the experiments. 

In the next section we briefly summarize the thermo- 
dynamic properties of the nodal p-wave superconduc- 
tor with A(k) oc sm(k x a) ± isui(k y a). In many re- 
spects the result* are very similar to the ones for d-waw 
superconductors^ and 3D /-wave superconductorslia. 
Then we proceed to consider the thermal conductivity 
in a planar magnetic field. The result for the 2D /-wave 
state with angular dependence cos(20)e ±I ^ is very similar 
to the one in d-wave superconductors discussed recently 
in Ref. Ill 



II. THERMODYNAMICS OF THE NODAL 
P-WAVE SUPERCONDUCTOR 

We consider the superconducting order parameter 
given by A(fc) = dj^[sm(k x a) ± ism(k y a)] with k x a — 
7r cos((f>) and k y a — 7rsin(</>) and the normalization sm — 
\/2sin(-^) = 1.125. This is the model proposed in 
Ref. ^ except that we have chosen the Fermi wavevector 
kpa = 7T in order to have a node in A(fc). The quasi- 
particle Green function in Nambu representation is given 
by 



G{k, lj) = (iuj - CfcP3 - A(fc)picri) 1 
where A(fc) = -^-[sm(k x a) ± ism(k v a)]. 



(1) 



0.5 




o Bonalde et al. 

2D f— wave 

nodal p-wave 

3D f-wave 



0.5 1 

VT C 

FIG. 2. The superfmid density p 3 (T) as a function of T/T c 
for the 2D /-wave states (solid line) and the nodal p-wave 
state (dashed Jine) together with the experimental data by 
Bonalde et al.u (circles) and the 3D /-wave state considered 
in Ref. |l| (dotted line). 



Then the quasi-particle density of states is given by 
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N(E)/N = Re( 
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where f 2 (<j)) — s^/(l — cos (V2ir cos <f>) cos(V27rsiri0)), 
y = E/A, and <> denotes an angular average. The 
density of states is calculated and shown in Fig. || to- 
gether with the one for the 2D /-wave case. In particular 
for E/A <C 1, the density of states increases linearly as 
N(E)/N ~ 0.7162EVA, while in the 2D /-wave case it 
varies like N(E)/Nq ~ E/A. Otherwise the two curves 
look very similar. Here, the gap equation 



A" 1 = </ 2 } 



2\-l 



f 2 E 
dE{ ■ I )tanh( — ) (3) 



has been solved numerically In particular we find 
A(0)/T c = 2.00, which has to be compared with 2.14 
in the 2D /-wave case. 
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FIG. 3. The density of states for the 2D /-wave state (solid 
line) and the nodal p-wave state (dashed line) . 

The entropy S is obtained from 

S = -4/ dEN(E)[f\nf+(l-f)\n(l-f)} (4) 



POO 

= 4 / dEN(E) [/3E(1 + e^y 1 + ln(l + e~ fiE )] 
Jo 

with / being the Fermi function and (3 = 1/T. Then the 
specific heat C S (T) is given by 

,dS(T) 



C S {T)=T- 



dT 



C s (T)/jT has been shown in Fig. |. 

Finally, the superfluid density p s (T) is given by 



Ps(T) = 1 



PA 



'dE^lsec^ m\ 
N I 2 J 



(•5) 



(6) 



which behaves almost linearly in T and is shown in Fig. 
^. We note that at low temperatures an expansion of 
p s (T) leads to p s (T) = 1 - 2 In 2 x 0.7162^ + • • •. 
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III. THERMAL CONDUCTIVITY TENSOR IN 
THE A-B PLANE 

As shown in earlier experiments on YBCO, the thermal 
conductivity tensor in a planat,iiiagnetic field is very sen- 
sitive to the nodal directionslijO and thus may be used 
to further discriminate between the states studied above. 
We shall consider the thermal conductivity tensor in the 
vortex state of the nodal p-wave and the 2D /-wave sepa- 
rately. We will show that only one of these states appears 
to be consistent with the angular independence observed 
recentlyEj'E-i 



B. 2D /-wave state 

As already mentioned we consider the two states 
sin(2</>)e ±l< ^ and cos(20)e ±Ic ^. The order parameter oc 
sin(20)e ± ^ has the same nodal structure as the nodal 
p-wave state studied in the last subsection and has been 
studied recently by Graf and BalatskyE]. Following the 
same procedure as above we find for the state sin(20)e ±l ^ 

k xx /k u = - <(1 + cos(20)) x) (x) = - V -^-F{9) (10) 
and 



A. Nodal p-wave state 

The necessary theoretical scheme, neglecting vortex 
core scattering, has been worked out during the past few 
yearalj'EJ'Ej. We just apply this method for the present 
case. In particular for < 1 in the superclean 

limit we obtain 
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where v and v' are the Fermi velocities within the a- 
b plane and perpendicular to it, respectively, and x = 
|v ■ q|/A denotes the Doppler shift due to the superflow 
around the vortex (see Ref. [TJ|) . n n — ^ r ^" is the normal 
state thermal conductivity. The function F(9) is given 
by 



F(0) = -Wl + sin 2 9E 




(8) 



with E being the complete elliptic integral of the second 
kind and 



1 x v 



/«n = 



(9) 



In the present situation there will be no off-diagonal com- 
ponent, because the heat current is parallel to the nodal 
direction. The angular dependence of n xx is given by the 
function F(6), which is shown in Fig. ^. Surprisingly, k xx 
has a broad maximum for 9 — n/2. Also, the anisotropy 
k xx (tt/2)/k xx (0) — 1.910 is quite strong. Therefore 
view of the thermal conductivity experimental datatJ 
we have to reject this possibility. 
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where F(9) has been shown in Fig. 
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FIG. 4. Angular variation of the functions F(9), I 2 (0), and 
I(0)L(0) as defined in Eqs. (|), Q, and @. The angular 
variation of the thermal conductivity Kxx for the nodal p-wave 
state, as given by the function F(6), is much stronger than 
for the cos(2^)e ±i</ ' /-wave state {I 2 (8)) due to the different 
position of the nodes in the gap function. 



Therefore, also the state sin(2(/))e ± ^ gives a rather 
large 9 dependence-j-which is inconsistent with the ex- 
istent experiments! 1 tF 4 !. 

Finally, let us consider the state cos(20)e ±l< ^, which 
has its nodes along the zone diagonal. As already noted, 
this state has the same thermodynamics as a d-wave su- 
perconductor. Eurther, the thermal conductivity tensor 
is now given byuJ 



and 



_ 2 vv'eH 2 
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where 
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and 



L(9) = 




(15) 



with s = sin(20). 

In Fig. H the angular dependences of the func- 
tions (1(0)) and I(9)L(9) are shown together with 
F(9). Thus, as in <i-wave superconductors, this state 
exhibits a fourfold symmetry in k xx . But the angu- 
lar dependence is ahe4it,10% and may be compatible 
with the experimentsliai-l Thus, we conclude that the 
cos(20)e ±Z1 ^ 2D /-wave state is the best candidate to de- 
scribe all of these experimental observations. 

As mentioned above, our analysis of the thermal con- 
ductivity tensor neglects vortex core scattering. At least 
in the high-T c compounds in a small magnetic field apd-al 
low temperatures this contribution can be neglectedtaO. 
In Sr2Ru04 the vortex core size is larger and at present it 
is unclear to what extend this contribution plays a role. 
We expect that the angular dependences shown in Fig. 
|] will be weakened both by vortex core scattering and 
finite temperatures, which will improve agreement with 
the experiments. 

As an additional check on the position of the nodes of 
the order parameter we propose a measurement of the 
transverse thermal conductivity n xy . As Eqs. (^|) and 
( |Tl| ) show, K xy vanishes, if the nodes lie along the a or b 
directions. However, we expect a finite transverse ther- 
mal conductivity n xy showing a sin(2#) variation for the 
cos(20)e ±,( ^ /-wave state having its nodes along the zone 
diagonal, as Eq. (O) shows. 



IV. CONCLUSIONS 

We compared one p-wave and two 2D /-wave super- 
conducting states with recent experimental data from 
purest crystals of Sr2Ru04. We find that within weak- 
coupling theory the two 2D /-wave states give the clos- 
est description of the thermodynamic data. Among 
these, the angular dependence of magnetotransport fa- 
vors the cos(20)e ±l( ^ /-wave state, since the other two 
states exhibit much stronger anisotropy than observed 
experimentally. Therefore, among the simplest states the 
cos(20)e ±l * /-wave state, having Bi ff x E u symmetry ap- 
pears to be the best candidate for superconductivity in 
Sr2Ru04. 
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